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I. INTRODUCTION 

The description of the hadron structure is related with 
our understanding of the non-perturbative properties of 
QCD. The essential part of the information about this 
structure is contained in the matrix elements of conserved 
quark operators. 

The electromagnetic current of a nucleon is 

J„ (P ,s';P, s) =u(P, 8 ) A M (q, P)u(P, s) 
= u(F', S ')(7 M Fi(<Z 2 )) + ^ia^q v F 2 (q 2 ))u(P, s), (1) 

where P,s,(P , s ) are the four-momentum and polariza- 
tion of the incoming (outgoing) nucleon and q = P — P 
is the momentum transfer. The quantity A M (g,P) is the 
nucleon-photon vertex. 

The similar expressions for the decomposition of ma- 
trix elements of energy momentum tensors describe the 
partition of angular momenta [l| and coupling of quarks 
and gluons to classical gravity @. The latter connec- 
tion manifests also relation to the equivalence principle 
for spin-gravity interactions [!, 0, [H, 0| and its possible 
validity separately for quarks and gluons 0, § • 

Two important combinations of the Dirac and Pauli 
form factors are the so-called Sachs form factors 
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with (t = —q 2 < 0). Their three-dimensional Fourier 
transform provides the electric-charge-density and the 
magnetic-current-density distribution [lCj| . Normaliza- 
tion requires G P E (0) = 1, G E (0) = corresponding to 
proton and neutron electric charges; Gm(0) = (Ge(0) + 
k) = fx defines the proton and neutron magnetic mo- 
ments. Here /i p = (1 + l.79)^j is the proton magnetic 
moment and k = F 2 (0) is the anomalous magnetic mo- 
ment: k v = 1.79 



Early experiments at modest t, based on the Rosen- 
bluth separation method, suggested the scaling behavior 
of both proton form factors and the neutron magnetic 
form factor approximately described by a dipole form 
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which leads to 



n 4M 2 -tu v 



4k p M 2 
(AM 2 - t] 



G L 



(4) 



(5) 



(0) 



with A 2 = 0.71 GeV 2 . These experiments were based on 
the Rosenbluth formula [ll| 
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where r = Q 2 /U1 2 and e = [1 + 2(1 + t) tan 2 (6» e /2)]- 1 
is a measure of the virtual photon polarization. 

Recently, better data have been obtained by using of 
the polarization method [H, [HI . Measuring both trans- 
verse and longitudinal components of the recoil proton 
polarization in the electron scattering plane, the data on 
the ratio 



r<p 



P t E + E 
Pi 2M P 
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were obtained. These data manifested a strong devia- 
tion from the scaling law and, consequently, disagreement 
with data obtained by the Rosenbluth technique. The 
results consist in an almost linear decrease of G P E /G P M . 
There were attempts to solve the problem by inclusion 
of additional radiative correction terms related to two- 
photon exchan ge a pproximations, for example [l4|. In re- 
cent works [TBI Il6l| the box amplitude is calculated when 
the intermediate state is a proton or the A-resonance. 
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The results of the numerical estimation show that the 
present calculation of radiative corrections can bring into 
better agreement the conflicting experimental results on 
proton electromagnetic form factors. Note, however, that 
the precise data of a Rosenbluth measurement of the pro- 
ton form factors at Q 2 4.10 GeV 2 [17( lie so high that 
they require very large corrections to move them down 
to meet the polarization data. 

The form factors are related to the first moments of the 
Generalized Parton distributions (GPDs) 0, OJ, OJ, HI • 
Generally, GPDs depend on the momentum transfer t, 
the average momentum fraction x — 0.5(xi + Xf) of the 
active quark, and the skewness parameter 2£ — Xf — Xi 
measures the longitudinal momentum transfer. One can 
choose th e sp ecial case £ = of the nonforward parton 
densities [2l| !F^(x\t) for which the emitted and reab- 
sorbed partons carry the same momentum fractions: 



n q (x,t) = H q (x,0,t) - H q (-x,0,t), 



(9) 



£ q (x,t) = E q (x,0,t) - E q (-x,0,t). (10) 
The form factors can be represented as moments 



F q {t) = dx H q (x,t), 
Jo 



F q {t)= [ dx£ q (x, 
Jo 



t), 



(11) 
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following from the sum rules P, [l9| . 

Non-forward parton densities also provide information 
about the distribution of the parton in impact parame- 
ter space [22| which is connected with i-dependence of 
GPDs. Now we cannot obtain this dependence from 
the first principles, but it must be obtained from the 
phenomenological description with GPDs of the nucleon 
electromagnetic form factors. 

Choosing a frame where the momentum transfer r is 
purely transverse r = rj_, the two-body contribution to 
the form factor can be written as I23| 



F«(i)= [ dz [ **(x,k ± + xr ± )*(x,k ± ) ^ 
Jo Jo 16tt- 
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As it was shown in [2l|, |24|, assuming the Gaussian 
ansatz 



V(x;k±) 
one can obtain 

F^(q 2 )-- 



exp[ 
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where (x) has the meaning of the two-body part of 
the quark density q(x). The scale A 2 characterizes the av- 
erage transverse momentum of the valence quarks in the 



nucleon and q 2 = —t is momentum transfer. This Gaus- 
sian ansatz was used in [2l[ for describing the form fac- 
tors of proton. However, this ansatz leads to a faster de- 
crease in Fi at large momentum transfer. So there arises 
an important question about i-dependence of GPDs. 

In [25| , it was proposed to use the factorized Regge-like 
picture 



H*(x,t) ~-^q(x) 



(16) 



which, however |24j, does not satisfy some conditions in 
the li ght- cone representation. 

In |26l. two parts of ^(x,k±) were introduced 
which are responsible for the interaction at small mo- 
mentum transfer ^(x, k±) so ^ t and for the interactions at 
large momentum transfer ^(x, k±) hard . These functions 
can reproduce the i-dependence of , for example, F\ but 
at the cost of the complicated mechanism and growth of 
F\ at hig her value of t . 

In [24[ , the Regge-like picture was used again but now 
without factorization 



H q {x,t) 



1 
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q(x). 



(17) 



This ansatz reproduces the basic properties of F\ and 
provides a good description of the ratio of F\ to F 2 
and, consequently, of the ratio of Ge to Gm- A sim- 
ilar approach was examined in [28| with GPDs form 
~ (x/go)~ a ^ t ^ 1 ~ x \ where a(t) is the nonlinear part of 
the Regge trajectory. 

Note that in [29|, [3(| it was shown that at large x — > 1 
and momentum transfer the behavior of GPDs requires 
a larger power of (1 — x) in the t-dependent exponent: 



H q {x, t) - exp[a (1 - x) n t] q(x). 



(18) 



with n > 2. It was noted that n = 2 naturally leads to 
the Drell-Yan-West duality between parton distributions 
at large x and the form factors. 

In other works (see e.g. [3l|, [HI, [33| ) the description of 
the ^-dependence of the GPDs was developed in a some 
more complicated picture using the polynomial forms: 



H q (x,t) ~q(x) Exp[f(x) t], 



(19) 



with 



f(x) = a'(l - x) 2 ln(-) + B q (l - x) 2 + A q x{l ~ a) (20) 



f{x) = a'(l - xfln(-) + B g (l - x) 3 + A q x(l - x) 2 (21) 

These parameterizations provide the good description 
of the data extracted from [34| The curve from Fig. 6 in 
[3ll | is reproduced in Fig.l. Note, that the parametriza- 
tion of (forward) parton distribution q(x) there also dif- 
fers from the other curves at Fig. 1. 
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II. NEW MOMENTUM TRANSFER 
DEPENDENCE OF GPDS 

Let us modify the original Gaussian ansatz in order to 
incorporate the observations of [2l[ and [3(| and choose 
the i-dependence of GPDs in the form 

H*(x,t) = q(x) exp[a + { ^-^- t}. (22) 

The value of the parameter m = 0.4 is fixed by the low 
t experimental data while the free parameters a± (a + 
- for 7i and a_ - for £ ) were chosen to reproduce the 
experimental data in the whole t region. Indeed, large t 
behavior corresponds to x ~ 1 in (10), (11), where the 
dependence on m is weak. 

The function q(x) was chosen at the same scale /x 2 = 1 
as in [H, which is based on the MRST2002 global fit 
In all our calculations we restrict ourselves, as in other 
quoted works, to the contributions of u and d quarks. 

Hence, we have 

u{x) = 0.262x-°- 69 (l - x) 3 50 (l + 3.83x a5 + 37.65z)(23) 
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FIG. 1: Proton Dirac form factor multiplied by t 2 (hard line 
- the present work, dot-dashed line - 12111 ; long-dashed line - 
[2^ ] ; dashed line - [27j ; dotted line - [3ll] ; the data for Ff are 
from [H| 



d(x) = 0.061aT°- 65 (l - x) 4 03 (l + 49.05a; - 5 + 8.65a;)(24) 

Following the standard representation, see for example 
(24| . we have for the Pauli form factor Fi 

£«(x,t) = £*{x) exp[a_ t}. (25) 

with 

e»(x) =^(l-xr u(x), (26) 

£ d (x) ^^(l-a;)^ d(x), (27) 

where K\ = 1.53 and K2 = 0.31 24]. 

According to the normalization of the Sachs form fac- 
tors, we have 

k u = 1.673, k d = -2.033, N u = 1.53, N d = 0.946 

The parameters a + = 1.1 and a_ were chosen to obtain 
two possible forms of the ratio of the Pauli and Dirac form 
factors. Below we consider variant (J) with a_ = 1.1 and 
variant (II) with a_ = 1.4 which correspond to the de- 
scription of the experimental data obtained by the polar- 
ization method and the Rosenbluth method, respectively. 

III. PROTON FORM FACTORS 

The proton Dirac form factor, multiplied by t 2 , is 
shown in Fig.l in comparison with other works and the 
experimental data. It can be seen that the complicated 
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FIG. 2: Ratio of the Pauli to Dirac proton form factor multi- 
plied by t (hard and dot-dashed lines correspond to variant (I) 
and (I I)) of the present work , dotted line - |46ll ; long-dashed 
line - [3) ; the data are from [47l. lij. |49|. |5(|. 



mechanism proposed in [26J leads to an increase in these 
values at higher t. One can see that our model repro- 
duces sufficiently well the behaviour of the experimental 
data at both high t and low t. 

The ratio of the Pauli to the Dirac proton form fac- 
tors multiplied by t is shown in Fig. 2. As it was men- 
tioned above, there are two different sets of experimental 
data. Here we are not going to discuss these two methods 
and different corrections to them. We just observe that 
our model describes the results of both the methods by 
changing the slope of £ only and present two respective 
variants of this slope a± . 

This is supported also by Fig. 3 where fi p G p E /G p M are 
shown in comparison with the experimental data. 
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FIG. 3: [i p G p E / G\ t (hard and dot-dashed lines correspond to 
variant (I) and (II)); the experimental data from [l7l. I48L |49|. 
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IV. NEUTRON FORM FACTORS 

Let us now calculate the neutron form factors using 
the model developed for the proton. The isotopic invari- 
ance can be used to relate the proton and neutron GPDs. 
Hence, we do not change any parameter and keep the 
same i-dependence of GPDs as in the case of proton 
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FIG. 4: G% (hard and dot-dashed lines correspond to variant 
(I) and (II)); the experimental data from [5g, |56|. |57|. 
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H n (x,t) = q+(x) n exp[2a+ (1 p f t], (28) 



where 



< \ 2 A 1 



(29) 



For the Pauli form factors of the neutron we corre- 
spondingly have 

£ n (x,t) = q.(x) n exp[2a_ t], (30) 

with 

q - {x)n = ~\w u {l - x)K1 u[x)+ \ih {l - xT2 d{xm 

We take the two values of the slope a_ as in the case of 
the proton form factors, corresponding to variant (J) and 
variant (//) below. 

First, let us calculate G%. The results are shown in 
Fig. 4. Evidently, the first variant is in better agreement 
with the experimental data. 

A more clear situation with our calculations of 
is shown in Fig. 5. In this case, it is obvious that the 
first variant much better describes the experimental data, 
especially at low momentum transfer. 

Finally, we present our calculations for the ratio of 
Ge/G 7 ^ for neutron in Fig. 6. 



CHARGE DENSITIES OF NEUTRON 



FIG. 5: GJf (hard and dot-dashed lines correspond to variant 
(I) and (II)); the experimental data from [581 ] . 




-t (GeV2) 

FIG. 6: GJ/Gm (hard and dot-dashed lines corresponds to 
variant (I) and (II)); experimental data from [HEI. |56| . 



motivated by recent discussion of the definition of charge 
density of the neutron at small impact parameters corre- 
sponding to the " center" of the neutron [3^, H3] • In , 
the charge density of the neutron is related to F™(t) and 
calculated using phenomcnological representation of the 
G%(t) and G n M (t). 



P(b) = 



^ e 9 J dx q(x, b) 



d 2 qF 1 (Q 2 = q 2 )e lirb 

" ^IjJqb) G E(q 2 )+rG M (g 2 ) 
2n 1 + t 



(32) 



Let us now discuss the neutron structure in the impact 
parameter (36l. l37l. |38| representation, we are particularly 



Jo being a cylindrical Bessel function. It differs essen- 
tially from the definition of the neutron charge distribu- 
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tion in the Breit frame related G%if). 

PG E (b) = J d 2 q[F 1 (q 2 )+TF 2 (q 2 )} 

= ^Mqb) G E (q 2 ). (33) 
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FIG. 7: Parton charge density of the neutron p n (b) corres- 
pound to Ge(P) 
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FIG. 8: Parton charge density of the neutron pn{b) corre- 
spond to Fi(b) 

Using our model of i-dependence of GPDs we may cal- 
culate both forms of the neutron charge distribution in 
the impact parameter representation and, moreover, de- 
termine separate contributions of u and d quarks. 

The charge distribution p n (b) corresponding to G% is 
shown in Fig. 7. It practically coincides with the result 
[4lj . The charge density of the neutron p n (b) correspond- 
ing to Fi is shown in Fig. 8. while the respective sepa- 
rate contributions of u and d-quarks are shown in Fig. 9. 
We can see that u-quarks have the large negative charge 
density in the centre of the neutron. 

VI. TRANSVERSE CHARGE DENSITIES 

Recently, the impact parameter distributions were gen- 
eralized [42| to the case of spin-flip magnetic form factor 
so that in addition to expression (|32f one has 

1 r°° n 2 

/4(b) = f$(b) + Sin{<t>)- / dq-l- Ji( ? &)F 2 (g a X34) 
2ir J 2M N 
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FIG. 9: w-quark (hard line) and (— d)-quark (dashed line) 
density of the neutron 



here tan((/>) = b x /b y . This transverse charge density for 
proton obtained in the framework of our model is shown 
in Figs 10 while the angular-dependent contribution (sec- 
ond term in (|3~i|) ) is shown in Fig. 11. 

The angular dependent part of the density is small and 
it is practically invisible in Fig. 10. However, removing 
the axially symmetric part of the density related to Pq (b) 
one can see that non-symmetric part is reaching the value 
about ±0.05 (see Fig. 11) at the distance of 0.3 fm from 
the center of the proton, comparable to the size of the 
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FIG. 10: Quark transverse charge density of the proton 
eq.(34): a)upper panel - back side view; b)low panel - view 
from above. 
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FIG. 11: Angular dependent contribution to quark transverse 
charge density of the proton p^(b) — PqQ)): a)upper panel - 
back side view; b)low panel - view from above. 

valence quark. 

The transverse densities of the neutron presented in 
Fig. 12 and Fig. 13. In this case the symmetric part of 
the transverse density p$ (b) is small (compare with Fig. 
8) . while the non-symmetric part is sufficiently large and 
concentrated around the center of the neutron. 



VII. GRAVITATIONAL FORM FACTORS 

Taking the matrix elements of energy-momentum ten- 
sor T^u instead of the electromagnetic current one can 
obtain the gravitational form factors of quarks which are 
related to the second, rather than the first moments of 
GPDs 

f dxxH q (x,A 2 ,0=A q (A 2 ); 

J dx x E g (x, A 2 , £) = B q (A 2 ). (35) 

For £ = one has 

/ dx xH q (x,t) — A q (t); / dx x£ q {x,t) = B q (t).(36) 
Jo Jo 

This representation combined with our model (we use 
here the first variant of parameters describing the ex- 
perimental data obtained by the polarization method) 
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FIG. 12: Quark transverse charge densities of the neutron 
po{b): a)upper panel - with the back side view point; b)low 
panel - with the upper view point. 
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FIG. 13: Quark transverse charge densities of the neutron 
pr{b) — po (b)- a)upper panel - with the back side view point; 
b)low panel - with the upper view point. 

allows to calculate the gravitational form factors of va- 
lence quarks and their contribution (being just their sum) 
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FIG. 14: Comparison of gravitational form factor A u+ d (hard 
line) and Proton Dirac form factor (dashed line) multiplied 
by t 2 , the data for F[ are from [45| 



to gravitational form factors of nucleon. Our result for 
A u +d{t) is shown in Fig. 14. Separate contributions of the 
u and d- quark distribution are shown in Fig. 15. At t = 
these contributions equal A u (t = 0) = 0.35 and Ad(t = 
0) = 0.14. The corresponding calculations for B q (t) are 
shown in Figs. 16. At t = these contributions equal 
B u (t = 0) = 0.22 and B d (t = 0) = -0.27. Hence their 
sum shows the sort of compensation supporting the con- 
jecture @, @]about validity of the Equivalence Principle 
separately for quarks and gluons: B u+ d{t = 0) = —0.05. 




FIG. 17: Difference in the forms of charge density iq and 
"matter" density (^4) 




-t (GeV2) 

FIG. 15: Contributions of the u-quark (hard line) and d- 
quark (dashed line) to the gravitational form factor A q 

Note that nonperturbative analysis within the frame- 
work of the lattice OCD indicates that the net quark 
contribution to the anomalous gravitomagnetic moment 
B u+ d{0)is close to zero 0,13. 

Let us compare the distribution of electric charge and 
matter (that is, gravitational charge) in the nucleon. For 
that purpose we generalize (|3^|) in a straightforward way 
and introduce the gravitomagnetic transverse density 
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FIG. 18: Densities of gravitational form factor A 



with a matter density 
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In Fig. 17 we compare this matter density with the charge 
density (c.f. Section 5) for proton. The plots for the 
angular dependent part of transverse density Px r {b) — 
Po r (b) are shown in fig. 18 and fig. 19. One can see that 
it is quite small and its maximal values are concentrated 
almost at the same distances from the center of nucleon 
as the transverse charge density of the proton. 
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FIG. 19: Transverse density p$ r (b) - p^ r {b) 

VIII. CONCLUSION 

We introduced a simple new form of the GPDs t- 
dependence based on the Gaussian ansatz correspond- 
ing to that of the wave function of the hadron. It satis- 
fies the conditions of the non-factorization, introduced by 
Radyushkin, and the Burkhardt condition on the power 
of (1 — x) n in the exponential form of the i-dependence. 
With this simple form we obtained a good description 
of the proton electromagnetic form factors. Using the 



isotopic invariance we obtained also a good descriptions 
of the neutron Sachs form factors without changing any 
parameters. We showed that both sets of the experimen- 
tal data (obtained by the Rosenbluth and polarization 
methods) on the electromagnetic form factors can be de- 
scribed by changing only the slope of the i-dependcncc 
of the spin-dependent GPD £(x,t). The comparison of 
these two variants with the neutron form factors gives 
preference to the one which describes the data obtained 
by the polarization method. 

Our calculations of the charge distribution of the neu- 
tron in the impact parameter form of F\ coincide with the 
calculation by Miller obtained from the phenomenologi- 
cal forms of G P E and G P M . They confirm that respective 
charge density is negative at small impact parameters. 
On the basis of our results we calculated the contribution 
of the u and d quarks to the gravitational form factor of 
the nucleons. The cancellation of these contributions at 
t = shows that the gravitomagnetic form factor is close 
to zero for separate contributions of gluons and quarks, 
which supports the conjecture of 0, H|. 
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